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ABSTRACT 

We present deep H I observations of the moderately inclined spiral galaxy, NGC 2997. 
The goal of these observations was to search for H I clouds in the vicinity of NGC 2997 
analogous to the high velocity clouds of the Milky Way and gain insight into their 
origins. We find evidence for the presence of a galactic fountain as well as the accretion 
of intragalactic material, however we do not identify any large clouds of H I far from 
the disk of the galaxy. NGC 2997 has a thick, lagging H I disk that is modeled with a 
vertical velocity gradient of 18 — 31 kms~^ kpc~^. Anomalous velocity H I clouds with 
masses of order 10^ Mq, which cannot be explained by galactic fountain models allow 
us to estimate a lower limit to the accretion of extragalactic gas of 1.2 Mq yr~^. The 
number and mass of these clouds have implications for cosmological simulations of large 
scale structure and the presence of dark matter halos. We have used values from the 
the literature to estimate a star formation rate of 5 ± 1 Mq yr~^ and to derive a new 
distance to NGC 2997 of 12.2 it 0.9 Mpc using published Tully-Fisher relations. 



Subject headings: galaxies: individual (NGC 2297) - galaxies: ISM - galaxies: kine- 
matics and dynamics - galaxies: structure 



- 2 - 



Introduction 



In the Milky Way Galaxy, 21 cm observations of neutral atomic h ydrogen (H I ) prov ide a 
unique probe of the structure and dynamics of the interstellar medium. iMuller et alj (|l963l ) and 
others have identified clumps of H I emission with motions that are inconsistent with Galactic 
rotation. Formally, these high velocity clouds (HVCs) are identified by either having a line of 
sight velocity, \vlsr\, greater than 90 kms~^ with respect to the local standard of rest or else 
differ fr om the rotational velocity of the disk by greater than 50 kms~^at a given Galactic lon- 
gitude (IWakker &: van WoerdenI 119971 ). Despite over 40 years of study, the origin of these HVCs 
remains a mystery. Our position within the Milky Way's disk provides a unique perspective of 
their spatial distribution. However, because the clouds do not obey Galactic rotation, distance 
measurement s rely on the serendipitous position of stars and st ellar spectroscopy to find distance 
brackets (e.g. iThom et al.ll2006l . 120071 : IWakker et al.l 120071 . 12008 ) or other coincidences, such as the 
associa tion of the Magellani c Stream with the Magellani c Clouds, or the interaction o f Smith's 
Cloud (|Lockman et ahlboosl ) and the Magellanic Stre am dMcClure-Griffiths et ahlbood) w ith the 



Galactic Disk. Other HVCs may be tidal debris (e.g. iLockman 



2003 



Put man et al 



200J). Such 



coincidences yield distances to individual HVCs of order 10 kpc, but are quite rare and tell us 
little about the population as a whole. Accurate distances to H I clouds allow us to calculate 
their masses and to map their positions about the Galaxy and thereby differentiate between several 
HVC formation scenarios, which in turn have implications for cosmology and structure formation, 
chemical enrichment, and halo dynamics. 

To explain the pervasive structure of Galactic HVCs across the sky, one must invoke multiple 
explanations. The Magellanic Stream is li kely the result of tidal inter actions w i th th e Large and 
Small Magellanic Clouds. For other HVCs. IShapiro &: Field (jl976l ) and lBregmanI (Il980l ) postulated 
that a galactic fountain mechanism could be responsible for raising hot gas into the Galact i c halo 
where it would cool, condense, and then rain back down on the disk. Meanwhile, lOortI ()1970l ) 
proposed that these H I clouds are primordial lef t-overs that are ju st no w being accreted by the 
Milky Way. This theory was revived recently by iBlitz et al.l (|l999l ) and iBraun &: BurtonI (jl999l ) 
who suggested that some types of HVCs were associat ed with low mass dark matter halos seen in 
models of ACDM galaxy f ormation (Klvpin et al.| 1999). HVCs may also be produced by a cooling . 



condensing hot h alo (e.g. iMaller fc Bullockl 12004 iFukugita Peebles 



2006 



Sommer-Larsen 



20061 . 



Peek et al.ll2008l ). These theories have different predictions for the mass distribution, the ver tical 
distance from the disk, and the metallicity for the clouds (see IWakker &: van WoerdenI 119971 . and 
references therein). But in either case, their presence would not be unique to our own Galaxy: in 
order to gain insight into the Galactic HVC population we observe other nearby spiral galaxies. 

As the closest large spi ral galaxy to the Mi lky Way, M 31 is an excellent test bed in which 
to search for HVC analogs. iThilker et al.l (|2004l ) find ^ 20 discrete objects within 50 kpc of the 
M 31 disk which range in H I mass between 10^ and 10'' Mq and can be characterized by a steep 
power law distribution. The authors suggest that these clouds may trace the dark matter halos 
around M 31 that are predicted from ACDM models. However, many of the clouds show velocity 
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continuity implying that they may be tidal remna nts associated with th e stellar stream of M 31, 
or part of the H I bridge between M 31 and M 33 (jWestmeier et al.ll2005l ). 



Observing edge-on galaxies reveals the vertical distribution of gas. The first, and subsequently, 



the m ost well studied galaxy for extraplanar H I is the edge-on spiral NGC 891 (e.g. ISwaters et al 
19971 ). The deepest publis hed observations prob e to H I masses of 9.3 x 10^ Mq and a column 



density of 1.6 x 10^^ cm~^ (jOosterloo et al.l 120071 ). These observations show that NGC 891 has a 
large H I halo that extends vertically out to 14 kpc from the thin H I disk and, in one quadrant, out 
to 22 kpc, contributing nearly 30% of the total H I mass of the galaxy. This halo is characterized 
by differential rotation, but it lags behind the rotation of the disk by 15 kms~^ kpc~^. Studies 
of NGC 891 and other nearly edge-on galaxies imply that there is a significant link between the 
amount of star formation in the disk (as measured by Ha or radio continuum), and the properties 
of the lagging halo, which is believed to be ev idence for gas circulation by a galactic fountain 
(|Swaters et all Il997l . iHeald et al.l boOBal fl l2007l v In general, they find that the scale height and 
the velocity gradient decrease with increasing star formation rate. Unfortunately, observations of 
edge-on galaxies are unable to tell us about either the position of clouds around the disk, or their 
vertical motion — the latter of which plays a part in identifying HVCs in the Milky Way. 

In face-on galaxi es H I at anomalous velocities and H I holes may be the result of gas ac- 
cretion fe.g. M 101: Ivan der Hulst Sancisil Il988l ) or massive star formation (e.g. NGC 628; 



Kamphuis Briggslll993 ). The mass and morphology of clouds are frequently used to discriminate 



between high velocity complexes that are likely to be in-falling extragalactic gas versus galactic 
fountain gas energized by supernovae. In the case of accretion, the gas tends to be located in the 
outer portions of the galaxy and is more massive. In the star formation case, the gas clouds tend 
to have lower inasses and are associated with regions of Ha or UV emission, or H I holes (e.g. 
Boomsma et al.ll2008l ). The largest H I holes in M 101 may be the result of recent collisions with 
extragalactic clouds; supernovae are ruled out because of the energy requirement to account for the 
amount of H I removed from the disk, and there is no evidence from UV or radio continuum data 



to support extraordinary star formation activity (Ivan der Hulst fc Sancisilll988l). NGC 6946 has an 
H I plume and large scale a symmetries probably related to tidal encounters (jKamphuis &: Sancisi 
19931 : iBoomsma et al.ll2008l ). as well as high-velocity gas and holes related to star formation. From 
the mass in the high velocity gas, iKamphuis &: Sancisil (119931 ) estimate a SN rate of 1 per 100 years 
to feed a galactic fountain, which is roughly consistent with the six SNe that have been observed 
in the galaxy this century. 

Inclined galaxies have the benefits of both edge-on and face-on observations. One can distin- 
guish the presence of extraplanar H I gas that lags behind the rotation of the disk and identify 
anomalous velocity gas with some sense of its location with respect to the disk while, at the same 
time, look for sites of active star formation that would be obscured in edge-on galaxies. There are a 
number of good examples of this. Position-velocity (hereafter p-v) slices of NGC 2403 show H I gas 
abov e and below the thin disk rotating 25 — 50 kms~^ slower ( Schaap et al.ll200d : iFraternali et al 



20021 ). analogous to the extra-planar gas in NGC 891. Furthermore, channel maps show counter- 
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rotating gas (iFraternali et al.ll200ll . 120021 '). In other examples, an as ymmetry in the distrib ution of 
diffuse gas in NGC 253 is evidence for merger or accretion activity (jBoomsma et al.ll2005l ). Mean- 
while Ha emission, froin star formation activity, is correlated with high velocity gas in NGC 6822 
(|de Blok &: Walteiil2006l ). H I ho les have been observ ed in NGC 4559 that are attributed to super- 
nova explosions or stellar winds ([Barbieri et al.ll2005l ). 



In order to better understand the origin of HVCs in the Milky Way and those seen around 
M 31, it is essential to observe a number of galaxies with similar properties that are well-situ ated for 



detai led study. NGC 2997 is the largest spiral galaxy i n the loose ga l axy g roup, LGG 180 ((Garcia 
19931 ) . This galaxy group was one of six observed by iPisano et al.l (|2007l ) in a Parkes survey of 



nearby Local Group analogs to search for intra-group HVC candidates. No HVC analogs were 
found in any of the groups down to a mass hmit < lO'^ Mq. We chose NGC 2997 for deep, high 
resolution observations because of its similarities to the Milky Way and M 31 in mass, luminosity, 
and star formation rate. It is a relatively isolated galaxy (no known companions within 110 kpc), 
it is moderately inclined with respect to the line of sight (~ 33°), and it has with a well behaved 
rotation curve. NGC 2997 was also chosen because it cannot be confused with Galactic or Local 
Group H I emission. Given a healthy star formation rate, it would not be surprising to find evidence 
for a lagging H I disk; its inclination allows us to investigate the position and kinematics of clouds 
with respect to the galactic disk in simple rotation. With these observations we identify extraplanar 
H I and HVC analogs associated with galactic fountains, galaxy formation, and tidal debris and 
use the characteristics of such HVCs to discriminate between their possible origins. 

The observations of NGC 2997 and data reduction are described in Section [2l The result- 
ing data are described in Section [3] including the results from modelling each H I component in 
NGC 2997. In Section [J] we consider possible origins of the extraplanar gas of the thick disk and 
the anomalous velocity H I clouds, and compare the H I in NGC 2997 to that of other galaxies. 
We also derive the star formation rate and distance of NGC 2997, and discuss how it affects the 
constraints on the origins of HVCs. We conclude in Section [5l 



2. Observations & Data Reduction 

In January 2006, we observed NGC 2997 with the Australia Telescope Compact Array (ATCA) 
for 28 hours in the EW682 configuration. These observations were in addition to archival ATCA 
data (from project C453) for NGC 2997 taken during 1995 and 1996, amounting to a total of 
59 hours in 6 different configura tions (Table [H) . The data were flagged and calibrated using the 



standard procedures in Miriad (jSault &: KileenI |2008| ). Prom 15-26 January 2007, we observed 



the galaxy for a total of 61 hours in 7-8 hour blocks with the Giant Metrewave Radio Telescope 
(GMRT). The GMRT is comprised of 30 antennas in a fixed "Y" -configuration with 14 antennas 
within 1 km and a maximum baseline of 25 km. The data were flagged and calibrated in AIPS 
at the GMRT using the standard procedures. Each night was individually imaged and inspected 
before being combined to create our final data cube. All observations were taken at a velocity 
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Table 1: Summary of ATCA Observations 





Configuration / 


Observing 


Date 


Max Baseline (m)^ 


Time (hrs) 


1995 Oct 27 


1.5D/1439 (4439) 


10 


1995 Nov 7 


6A/2923 (5939) 


14 


1995 Dec 15 


6C/2786 (6000) 


13 


1996 Jan 21 


750C/750 (5020) 


12 


1996 Jan 23 


750B/765 (4500) 


11 


1996 Feb 22 


1.5C/1485 (4500) 


9 


2006 Jan 7 


EW352/352 (4439) 


8.5 


2006 Jan 8 


EW352/352 (4439) 


8.5 


2006 Jan 9 


EW352/352 (4439) 


11 



^Mciximum separation between antennas, excluding (and including) the 6 km antenna. 

resolution of 6.596 kms~^. 

To combine the data, the reduced GMRT observations were written out of AIPS and read into 
Miriad 0. The visibility data were combined in Miriad with the task UVAVER, and cleaned using 
the task MOSSDI. We imaged 120 channels centered on the galaxy using a robustness of 2 and 
down-weighted the longest baselines from the GM RT. We derive a dist ance to NGC 2997 of 12.2 



Mpc using published 2MASS Tully-Fisher relations ([Masters et al.ll2008l ) that we adopt throughout 
the paper. Discussion of this derivation is left to the Section 14.31 in order to not take away from 
the focus of this paper on the H I properties of the galaxy. At this distance, our observations have 
a spatial resolution of 14" = 0.8 kpc. Altogether, our final cube has a rms noise of o" = 0.404 
mJy kms""*^ beam~^ channel^-*^. At 5a we achieve a mass sensitivity of 4.7 x 10^ Mq channel""^ 
and a column density of 9.6 x 10^^ cm~^ beam~^ channel"^. 



3. Results 

We present our data cube in the form of channel maps in Figures [T][3j These channel maps 
have been clipped at the 3a level to emphasize the galactic emission. The figures illustrate the 
general rotation of NGC 2997, but also show evidence of a thick disk and additional H I that is 



^In Miriad, the imported image headers were incomplete: we had to set the rest frequency to 
1.420405752 GHz and, for the purposes of estimating the rms noise before imaging, we added the cor- 
rect values for the system temperature (76 K) and system gain (0.22 Jy/K/antenna x 30 antennas). 
^http: //www.gmrt .ncra. tif r .res . in/ginrt_hpage/Users/doc/ spec_mod.pdf ^ 



- 6 - 



moving at velocities inconsistent with a rotating thin disk. 

In Figured] we present moment maps created from the data cubes: the optical DSS image with 
H I contours, the total H I intensity map, the velocity field, and the velocity dispersion (moments 
0-2). The moments are calculated independently for each pixel within a specified clip range — in 
our case, it is calculated from pixels containing emission above 5a of our rms sensitivity in order 
to avoid large outliers, particularly on the outskirts of the galaxy where the noise is increasing 
with distance from the pointing center. In the optical image, it becomes evident that the H I gas is 
extended to the northwest with respect to the optical disk. However, the intensity weighted velocity 
map looks fairly well behaved; there is no direct evidence that the galaxy has been disturbed by a 
strong encounter. Finally, the velocity dispersion map shows that the greatest dispersion occurs in 
the center likely due to beam smearing effects and due to intense star formation as suggested by 
strong radio continuum emission at the center of the galaxy. 

Figure [5] shows the integrated H I profile over the area of the entire galaxy. Summing the 
emission in this profile yields a total fiux for NG C 2997 of 226.6 Jy kms~^. This compares to 



162.3 Jy kms ^ from HIPASS ( Mever et al.ll2004l ). However, the latter value is likely to be low 



since the HIPASS measurement assumed the galaxy was unresolved. At a distance of 12.2 Mpc, 
our measured fiux corresponds to a total H I mass of 7.96 X 10^ Mq. 

We are interested in understanding the nature of the H I in NGC 2997 and identifying gas 
clouds that are analogous to the HVCs observed in the Milky Way. In decomposing the galaxy, 
we find there are three components that contribute to its total H I content. A thin disk dominates 
the total H I emission and can be modeled by a Gaussian profile with a velocity dispersion of 
11 lb 3 kms~^. We modeled and fit a rotation curve to this thin disk, however, there is additional 
gas that appears as a wing on the side of the thin disk's velocity profile towards the systemic 
velocity. This "thick disk" contributes approximately 16 — 17% of the total H I mass and has a 
velocity dispersion of ~ 16 kms^^. While this second component shows regular rotation, it lags 
behind the thin disk at a relatively constant rate. The third component consists of H I emission 
that cannot be explained in terms of either disk component, and some of which is counter-rotating. 
These three components stand out in Figure [H 

Table [2] summarizes the properties we derive for the galaxy. We ca l culate a new distance to 



NGC 2997 based on the 2MASS Tully-Fisher relations from lMasters et al.l (|2008l ). The radio center, 
systemic velocity, inclination, position angle, and the maximum rotational velocity come from our 
rotation curve fit to the velocity field. The integrated flux and Mhi come directly from the data 
cube. Mdyn is calculated from the maximum rotational velocity. The velocity gradient, scale height 
and cyclindrical height come from our best fit model of the extraplanar H I. The star formation 
rate is estimated from values in the literature. Sections I3.m3.3l and Discussion describe in detail 
how we arrive at these values. 
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Table 2: Properties of NGC 2997 



Property Value 



Hubble Type SAc^ 

Optical Center (J2000) 09^45'"38.8'^ -31°11'27.9"" 

Optical Extent 8.9'x 6.8"" 

Distance 12.2 Mpc"" 

Radio Center (J2000) 09''45™35.9" -31°11'27.8" 

Radio Extent 15.6' x 11.1'= 

Systemic Velocity 1089 it 5 kms^^ 

Inclination, i 32.3° ± 0.3 

Postion Angle 107.1° ± 0.6 

Integrated Flux Density 226.6 Jy kms~^ 

H I Mass, Mhi 7.96 x 10^ Mq^ 

Maximum Vrot 226 ± 2 kms^^ 

Dynamical Mass, Mdyn 2.1 x 10^^ M© 

Velocity gradient, ^ 18 — 31 kms""*^ kpc~^ 

Scale Height, z 0.7 ± 0.2 kpc 

Cyclindrical Height, z^yi 0.5 ±0.1 kpc 

Star Formation Rate 5 it 1 Mq yr~^ 



Note. — Unless otherwise stated, the values in Table [2] are derived from this work. 
"Optical position and extent obtained from the NASA/IPAC Extragalactic Database (NED). 
''Distance and H I mass are derived from the 2MASS TuUy Fisher relations of Masters et al. 1 2008h . 
'^Measured from 3a contours along major and minor axes. 



3.1. Modeling the Thin Disk 



As illustrated by Figure [U the H I profile towards most positions in NGC 2997 is dominated 
by a narrow velocity component with a dispersion of 11 it 3 kms^^. To m odel it, we clipped the 
data cube at 10% of the galaxy's peak emission (in a manner similar to iFraternali et al.l |2002| ) 
and fit a Gaussian to this upper part of the velocity profiles. By clipping the emission at this 
level we prevent the subsequent analysis of NGC 2997's thin disk from being biased by broad, low 
column dens ity H I ve l ocity components. We fit a rotation curve to the thin disk following the 
procedure of iBegemanI (|l989l ). We assume that the position angle and inclination of the galaxy are 
constant for all radii (we assume an ideal, flat disk), and from this we derive a systemic velocity 
of 1089 lb 5 kms~^, a position angle of 107.1° it 0.6 and an inclination of 32.3° it 0.3. The errors 
are the average errors found for rings between 5 — 30 kpc, however. Figure [7| shows that there is a 
gradient in the position angle, which has not been taken into account. 
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The top panel of Figure [7] shows our best fit rotation curve. It rises out to ~ 12 kpc, then 
slowly declines to greater than 25 kpc. We trust this rotation curve to ~ 28 kpc, beyond which the 
lopsidedness of the H I emission produces large errors in the fit. In addition, in order to understand 
the quality of our assumption that the disk is flat and ideal, we re-fit the rotation curve allowing the 
inclination and the position angle to vary while keeping the rotational velocity fixed to the value 
found for each ring in the tilted ring model. The middle and bottom panels of Figure [7] show how 
the inclination and position angle, respectively, vary with radius. In these instances the rotational 
velocity for each ring is fixed to the value plotted in the top panel. 

There is evidence that NGC 2997 deviates from the ideal, flat disk. The bottom panel of Figure 
[7] suggests that we were unable to fit a self-consistent rotation curve with a uniform position angle 
because it declines steadily with radius. In the tilted-ring model, inclination and rotation velocity 
are coupled through a Vrot sin i term. This makes for a degeneracy between a flat inclination and 
falling rotation curve that is particularly strong for low inclination galaxies. The "S" like appearance 
of contours near the minor axis (Figure H]) indicate that NGC 2997 may have a slightly warped H I 
disk. 



3.2. Modeling the Lagging Thick Disk 



Our deep observations of NGC 2997 show that the H I gas is not confined to the thin disk. 
In p-v slices the contours of H I emission show a shallow gradient towards the systemic velocity: a 
feature which has been id entified in a small number of other galaxi es and has been referred to as a 
"beard" (e.g. NGC 2403: ISchaap et al.l 120001 : iFraternah et al.ll2002l ). After modeling the dominant 
H I component, we subtract the thin disk from the original data cube and find that 1.4 x 10^ Mq 
of residual H I emission remains — about 18% of the total Mhi of the galaxy. 

This second component can be fit by a Gaussian with a higher velocity dispersion than the 
thin disk (~ 16 k ms~^) which reflects th e thermal motions of the gas convolved with a vertical 
velocity gradient (jde Blok k. Waited 120061 ) . The peak velocity of the fitted Gaussian resides closer 
to the systemic velocity of the galaxy than gas in the thin disk at the same radius. In other words, 
it appears to follow galactic rotation, but at slower rotational velocities. We fit a second rotation 
curve to the residual data by assuming that this thick disk has the same large-scale morphology as 
the thin disk, with the same systemic velocity (1089 kms~^), position angle (107°) and inclination 
(32°). From Figure [8] we see that peaks in the two rotation curves roughly coincide with one 
another, implying that they are related, however the thick disk consistently lags behind that of the 
thin disk by ^ 60 k ms~^. A similar lagging disk has been observed in NGC 891 (25 — 2 00 kms~^; 



Swaters et al.lll997l ) and NGC 2403 (25 - 50 kms~^: ISchaao et al.ll200d : IFraternah et al.ll2002l ') 



To understand this relation we also attempted to derive the thick disk rotation curve from 
the residual data without assuming it has the same properties as the thin disk. While we success- 
fully fit a position angle (~ 25°), we were unable to find a consistent inclination. In NGC 2403, 
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Fraternali et al.l (|2002l ) show that a rotated position angle, such as this can be explained by an 
overall inflow of the anomalous gas. However, as we find no other evidence for radial gas motion 
when we include it in our models of NGC 2997, we disregard it. 

We took this kinematic modeling a step further by adding a velocity gradient to the rotation 
curve as a functi on of height a bove t he dis k mid-plane. W e foltow ed similar procedures to those 
demonstrated bv iHeald et aD (|2006al l and lOosterloo et al.l ^20oi ). The GIPSY task GALMOD 
allows one to create a model data cube from user supplied parameters: the rotational velocity, 
the scale height and the surface density. A modified version of this task (Heald, 2008, private 
communication) includes two additional parameters to specify (1) the rate of change in velocity, 
dv/dz, with height, z, and (2) a cylindrical height above the disk at which the gradient begins, z^yi, 
such that the rotation of the gas is described by the following equations: 



Vrot{R,z) 



(1) 



VrotiR,0) Z < Zcyl, 1 

VrotiR, 0) - {z - Zcyl)^ Z > Zcyl J 

We examined p-v slices along the galaxy's major axis to compare the models created from the 
modified GALMOD task with our data (Figure [9|) . 

The cylindrical height was the easiest parameter to fit. It is defined as the distance from the 
mid-plane at which the velocity gradient begins. Setting its value has the eff'ect of imposing the 
velocity at which the peak of the emission occurs. For example, with a smaller cylindrical height 
the velocity gradient begins closer to the mid-plane and, for an inclined galaxy, the resulting sum 
of emission projected along the line of sight appears to shift towards lower velocities. A cylindrical 
height of O.S i liO.l k pc fits our data best. This is within the range of values derived for other galaxies: 
Heald et all pOOTI ) report values ra nging from 0.4 kpc in N GC 4302 to 1.2 kpc in NGC 5775 using 



(12003) report a value of 1.3 kpc for NGC 891 



Ha kinematics of ionized gas; while lOosterloo et al 
in H I. 

The value of the scale height predominately infiuences the lowest level H I contours of the p-v 
diagram. A large scale height places more emitting material at greater distances from the disk 
and, therefore, at lower velocities with respect to the rotational velocity at the mid-plane so the 
H I contours become shallower toward the systemic velocity with increasing scale height. The best 
fit for the scale height is 0.7 it 0.2 kpc. Our errors are based on the resolution of parameter space 
sampled and the degree to which the models differed significantly by eye. 

The model is least sensitive to the velocity gradient, probably due to NGC 2997's relatively low 
inclination and the low column density of the gas in the lagging disk. Varying this value changes 
the gradient of the H I contours in the model, however it does so less dramatically than varying 
the scale height. This, combined with the dumpiness of the H I due to spiral arms and other 
features in the disk, make it difficult to estimate a good fit. We estimate the velocity gradient to 
be between 18 — 31 kms~ ^ kpc~^. Again, these values fall in the range reported for NGC 4302 



NGC 891 (15 kms-i kpc"^: lOosterloo et al.ll2007l ^. 



(22.8-59.2 kms-i kpc~^: lHeald et al.ll2007l). NGC 5775 (8 kms'^ kpc'^ |HeaId_etalJl2006a|), and 
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3.3. Anomalous H I 

In spite of this extensive modeling, there is a significant amount of remaining H I that is 
spatiahy or kinematicahy distinct from the disk gas. In order to identify this faint emission we 
Hanning smoothed in velocity and spatially smoothed the data to twice the original beam size. The 
resulting cube has a resolution of 28"x22". The rms noise improves to cr = 0.15 mJy beam~^, or 
a 5(7 mass sensitivity of 1.7 x 10^ Mq channel"^ and column density of 8.9 x 10^^ cm~^ beam"^ 
channel"^. We visually searched this cube and found several examples of anomalous H I features 
whose basic properties are presented in Table their projected radial distance from the center of 
the galaxy, their velocity offset from the H I disk, and their H I mass as well as an estimate for 
their contribution to the mass accretion rate. This is not an exhaustive list, but merely illustrates 
some of the most prominant examples of HVC analogs associated with NGC 2997. Figure [TO] shows 
where we have taken slices through the smoothed data cube to present these features in Figures 
dUA-D. 

Of these anomalous features, we find 5 examples of H I clouds that are completely distinct from 
the H I disk. They were identified in p-v slices as having velocity offsets of greater than 50 kms~^ 
from H I gas in the disk at the same position (Figures 1101 fc fTT]) . In three of the five examples, the 
H I gas is counter-rotating with respect to the H I disk (Slices A & C), and in one instance we find 
an H I cloud whose velocity component along the line of sight is actually faster than the rotation 
of the disk (Slice B). The clouds we identify reside within 20 kpc from the center of the galaxy, 
which places them within the extent of the projected H I disk, but outside the optical disk as it is 
seen in DSS im ages (Figure B l ), and they all have H 1 masses of approximately 10^ Mq. They were 



not detected by lPisano et al.l (j2007l ) because of the low angular resolution of their Parkes (14') and 
ATCA (~ 1') data. There are significant indications of a larger population of H I clouds around 
NGC 2997, however we are less confident in their identification because they are closer to our mass 
sensitivity limit and because of residual phase calibration errors in the data cube. 

In addition to the H I clouds, there are two other prominent features of anomalous H I we see 
in velocity channels: a "spur" of counter-rotating gas near the center of the galaxy and a "tail" 
of emission lagging in rotation, extending to the west from the north half of the galaxy. Unlike 
the H I clouds which are distinct, both of these features are connected to the main disk of the 
galaxy (making it difficult to estimate their mass). The "spur" is attached north of the center, 
and arcs to the east and south extending greater than 6 kpc east from the center of the galaxy 
(f = 1016.1 — 1082.0 kms^^ in Figures [U & [2]). In p-v slices the same feature looks like a shell 
of emission that is connected to the disk (Figure [TTB. C & D). We estimate the H I mass in the 
spur to be ~ 2.6 x lO"^ Mq. The "tail" contributes to the asymmetry total intensity map (Figure 
HI). Although the contours of the intensity weighted velocity map (Figure H]) are well behaved, 
individual velocity channels, v = 943.5 — 983.1 kms~^ (Figured]), show emission on the north side 
of the galaxy. The coherence of this emission can be seen in Slices A & B of Figure [11) the feature 
looks like a bifurcation of the H I disk in p-v space. Its velocity component along the line of sight 
is actually moving faster than the rotation of the disk. We calculate the mass of the tail to be 
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3.1 X 10^ M(T,. 



Table 3: Properties of Anomalous H I Clouds in NGC 2997 





Distance 


Velocity*^ 


TF Mass"" 


Accretion Rate'^ 




Obj. ID 


(kpc) 


(km s^"*^) 


(Mo) 




{Mq yr"^) Notes 


Cloud 1 


17.7 


120 


1.9 ±0.2 X 10"^ 


0.13 


Fig. dUA 


Cloud 2 


13.0 


-105 


7.6 ±0.1 X 10^ 


0.06 


Fig. 


Cloud 3 


18.9 


150 


8.6 ±0.4 X 10^ 


0.07 


Fig. duc 


Cloud 4 


14.2 


80 


8.9 ±0.8 X 10^ 


0.05 


Fig. ED 


Cloud 5 


10.4 


200 


4.9 ±0.4 X 10^ 


0.10 


Fig. dUA 


H I Spur 


5.9 


150 


2.6 ±0.1 X 10"^ 


0.70 


Fig. dlB, C & D 


H I Tail 


27.8 


-60 


3.1 ±0.3 X 10"^ 


0.07 


Fig. iniA. & B 



"Negative indicates a line of sight velocity faster than the rotational velocity of the disk at that position. 
''The accretion rate is estimated from Mhi cloud x 'Voffset/r where Voffset is the velocity difference between the cloud 
and gas in the thin H I disk, and r is the projected distance of the cloud from the center of the galaxy. 



4. Discussion 

In addition to the thin disk of NGC 2997, seen in most spiral galaxies, we have identified two 
additional H I components that are not as ubiquitous: a lagging, thick disk and anamolous H I 
clouds. The thick disk has a mass of 1.4 x 10^ Mq, a broader velocity dispersion, ~ 16 kms~^, and 
lags behind the rotational velocity of the thin disk by ~ 60 kms^^. The anomalous H I clouds are 
spatially or kinematically distinct from both H I disk components. The clouds have individual H I 
masses between 1O^'^~'''^M0, are offset from the disk by 80 — 200 kms~^ and, in some instances, 
are counter-rotating. In this section we compare our results with the distribution of H I observed 
around other galaxies, we speculate on the possible origins of NGC 2297's HVC analogs, we discuss 
our results in the context of the Local Group HVCs, and we examine the effect of our distance 
estimate on our conclusions. 



4.1. Extraplanar H I 

NGC 2997 is one of the lowest inclination galaxies in which a lagging disk has been modeled. 
It is impossible to do this same decomposition for a face-on galaxy because the measured velocities 
contain no information on the rotation of the gas, but moderate to high inclination galaxies are 
well suited to this analysis. The velocity gradient we estimate for NGC 2997, 18 — 31 kms~^ kpc~^, 
appears to be typical of spiral galaxies for which deep observations reveal extended extraplanar gas. 
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Our findings are consistent with that seen in NGC 2403 where iFraternali et alj (|200lh calculate 
the thick disk to hav e a mean rotat i on vel ocity 25 — 50 kms~^ lower than that of the thin disk, 
in NGC 4459 where iBarbieri et al.l (|2005l ) find a thick disk that lags by 12 — 25 kms~^, and 
in NGC 891, where t he velocity gradient r anges from 14 to 43 kms~^ kpc~^ depending on the 
location in the galaxy ([Oosterloo et al.ll2007l ). NGC 6946 also has a similar lagging halo of cold gas 
(jBoomsma et al.ll2008l ). 



This range of values for the velocity grad i ent in spi ral ga laxies suggests that one might correlate 
it with star formation activity. iHeald et alJ (|2006al lbl. 120071 ) have done a complementary study of 
three high inchnation galaxies, NGC 5775, NGC 891, and NGC 4302, in which they find that the 
magnitude of the velocity gradient, measured in ionized gas (Ha), decreases with star formation 
activity, as measured in far infrared luminosity. For NGC 2 997, we calculate a s tar forma tion rate 



of 5±1 M0 yr-i from 60 and 100 urn IRA S fluxes_(following iKewlev et al.ll2002l : lBellll2003l ) and 1.4 



GHz fluxes (following ICondon et al 



2002 



Belli I2OO3I ). However, it is unclear how weh NGC 2997 



supports such a correlation: its values for the velocity gradient a nd scale height are close to the 
mean of those galaxies, while the star formation rate is a bit high. iHeald et al.l (120071 ) also suggest 
that the product of the velocity gradient and the scale height tends towards 20 km s~^. We estimate 
dv/dhz = 12 — 21 kms~^, which is consistent with their suggestion. 

From a theor etical standpoin t, the detection of a thick, lagging disk is consistent with the 
idea developed by iBregmanI (Il980l ) in which the Galactic HVCs are clouds condensing out of a hot 
Galactic corona fed by supernovae in the disk. In his galactic fountain model, gas blown out of the 
disk by star formation moves radially outward and slows down to conserve angular momentum as it 
moves vertically away from the mid-plane. If the observed lag is only due to conservation of angular 
momentum then, in the case of NGC 2997, gas that lags by 60 kms~^ behind a disk rotating at 
200 kms~^ originated 70% closer to the center of the galaxy. However, if hydrodynamic effects 
are important, the radial motion is less dramatic. Cooling gas condenses out of the corona and 
moves inward again. E vidence for radial inflow of 1 — 20 kms~ ^ has been observed in NGC 2403 



(jFraternali et al.l 120011 ) and perhaps in NGC 4559 (IBarbieri et al.ll2005l ). We do not find evidence 
for radial inflow of H I in NGC 2997, but such a signature may be ambiguous due to the low 
inclination of the galaxy. Radially infalling clouds would also be expected if these clouds originated 
in a cooling hot galactic ha lo that is at the virial temperature of the dark matter halo as proposed 
by iMaller &: Bullockl (|2004l ) , however holes in the H I disk may be evidence that supernovae are 
responsible for circulating this gas. 

In a number of cases where we see anomalous H I clouds in p-v slices, we also see depressions 
in the H I column density in the disk (Slices B &: D). These H I "holes" occur near the outskirts 
of the H I disk, but outside of the optical disk, and we hesitate at calling them holes because they 
frequently occur between two H I spiral arms. These may be sites where H I has been blown out of 
the disk from supernovae, however one cannot rule out that these may also be the result of a dwarf 
galaxy that has impac ted the disk. For example, this may be the case for large holes in NGC 6946 
(jBoomsma et al.ll2008l ) and M 101 (jvan der Hulst &: Sancisilll988l ). In any case, it is interesting to 
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note that all of the anomalous H I we identify is coincident with the H I disk and we do not find 
any large H I clouds far from it. 

The large clouds of counter-rotating H I gas and the "spur" cannot be explained by a galactic 
fountain mechanism. While fountain gas lags behind the main disk in rotation, it will conserve 
angular momentum and so can not be counter-rotating. Additionally, to eject clouds of lO''' Mq 
from the disk at the observed velocity offsets requires an spatially and temporally coherent injection 
of 10^^ ergs of energy. This corresponds to of order 10^ supernovae, assuming 10% efficiency. Despite 
the star formation rate of NGC 2997, this appears to be implausibly high. Based on this energy 
argument, and their position with respect to the optical disk, we conclude that these H I features 



are m ore likely material being accreted by NGC 2997, either out of a hot galactic halo (|Peek et al 



20081 ). or are the re mnants of a t idally disrupted gas rich dwarf galaxy. The maximum individual 



cloud mass from the iPeek et al.l (|2008l ) simulations is 5 x 10 Mq, however they do not rule out 
the possibility of cloud complexes having higher masses. We estimate a lower limit for the total 
accretion rate of 1.2 Mq yr~^, which may help sustain the vigorous star formation rate, especially 
if it is augmented by gas from the lagging disk which one would expect to eventually fall back 
towards the mid-plane. 

Features such as the asymmetry in the extended northwest spiral arm are usually interpreted 
as evidence for past interactions between galaxies, and there are two possible culprits. The nearest 
LOG 180 group galaxy is a dwarf galaxy, ESO 434-G041, which is ~ 115 kpc to the southeast. 
It is resolved in our ATCA observations, but lies outside the primary beam of the GMRT. NED 
lists another dwarf galaxy, ESO 434-G030, ~ 90 kpc to the southwest, but we do not see this 
galaxy in our data. Assuming that the velocity differences in the plane of the sky are similar to 
the line-of-sight velocity differences between these galaxies and NGC 2997, any encounter would 
have taken place ~ 1 Gyr ago. While there are no apparent connections between these galaxies 
and NGC 2997 in H I, nor any optical signatures of disturbances in the dwarf galaxies, a past tidal 
event or interaction seems to be the most likely explanation for pulling such a large amount of H I 
from the disk in such a coherent manner, causing the asymmetry. On the other hand, we do not 
believe that interactions are the most probable origin for the anamolous H I clouds. Similarly, the 
"tail" feature, on the northwest side of the galaxy, seems best explained by accretion-perhaps a 
dwarf galaxy that has been tidally disrupted as it passed near NGC 2997. 



4.2. Comparison with M 31 Clouds and Galactic HVCs 

We have identified the most massive anomalous H I clouds associated with NGC 2297. These 
five clouds have masses of order 10 Mq. The spectrum of cloud masses extends to lower masses 
but, in many cases, it is difficult to separate them as distinct from emission in the thick disk. 
Nonetheless, we can still compare the clouds we have identified with those reported in the Local 
Group and around other galaxies. 
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The H I mass distribution of clouds in M 31 given by dM(M/f7) oc Mhi~'^(^M}{i, allows us to 
estimate the total number of similar clouds we should expect to see around NGC 2997. Most clouds 
have masses in the range (0.15-1.3) X 10^ Mc:, with only Davies' Cloud having Mm = 10^ Mq. For 
sources smaller than our beamsize, seen by iThilker et al.l (j2004l ) we would expect to detect about 
10 of the M 31 HVC analogs at the 5cj level (and about 4 above the lOo" level); for those clouds 
that we could resolve, we lack the column density sensitivity to detect them. However, we already 
know from the observations of lPisano et al.l (j2007l ) that there are no low column density H I clouds 
res iding off the disk of NG C 2997. The HVCs discovered around M 31 in higher resolution images 
by Westmeier et al. (j2005l ) have an average size of 1 kpc (equivalent to the size of our beam at the 
distance of NGC 2997) and are below our detection limits. 

With recently determined d istant brackets to Mi lky Way HVCs, we can put our H I clouds into 
a local context. For Complex C. lwakker et al.l (120071 ) estimate Mhi to be (0.7-6) x 10*^ Mq. The 
H I masses for four ot her clouds range from (0.09-0.41) x 10"^ Mq (Cloud gl) to (0.68-1.60) x 10^ 
Mq (Complex GCP; Iwakker et alJbood ). Therefore, at the distance of NGC 2997, we are only 
sensitive to the most massive Milky Way HVCs for which distance measurements exist. However, 
we also see more massive features (~ 10^ ^o) for which their are few counterparts around the Milky 
Way and M 31 . Nonetheless, our discoveries are consistent with anomalous H I clouds observed in 
other galaxies (jSancisi et al.ll2008l . and references therein). 



4.3. Implications of the Distance Estimate to NGC 2997 

The distribution of H I clouds we find around NGC 2997 can be used to place limits on 
cosmological and large scale structure simulations for the number of small dark matter halos around 
massive galaxies. The H I masses we calculate depend on the square of the assumed distance to the 
galaxy. However, local large scale structure results in an inaccurate distance estimate to NGC 2997 
because its peculiar velocity is not dominated by H ubble flow. Unfortunately, an attempt to correct 



for this, the velocity flow model of iMastera (|2005l ). is sparsely populated by galaxies in the volume 



around NGC 2997, so we choose to calculate the distance to NGC 2997 using the Tully-Fisher 
relation (hereafter TF). 

The TF relations based on near infrared observations may be the most reliable because they are 
less impacted by dust than optical magnitudes. We retrieved the 2MASS J, H, and K-band magni- 



tudes from NED and used the bias corrected TF parameters speciflc to Sc galaxies (jMasters et al 



20081 . Table 2): 



M( J) - 5 log = -21.017 - 9.228(log W - 2.5) 
M{H) - 5 log = -21.847 - 9.165(log W - 2.5) 
M{K) - 5 log /i = -22.039 - 10.092(log W - 2.5) 



To estimate the internal extinction, we interpolated between I and K bands to get values for J and 
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comes from Parkes observations ( Pisano et al 



H using the formulae in lCourteau et alj (12007). Th e FWHM of the H I profile {W = 445 kms~ 



20071 . Pisano et al. 2009, in preparation), corrected 



for instrur uental effects, t urbul ent motions, cosmological broadening, and inclination as per the 
formula in [Masters et al.l (|2008l ). From these 2MASS TF relations we derive a mean distance of 
12.2 ±0.9 Mpc. 

This new distance estimate places NGC 2997 about 10% closer to the Milky Way than the 
estimates found on NED, corrected for peculiar velocit y due to local density enhancements, and 
18% closer than the velocity flow model of lMastersI (|2005l ). This places somewhat tighter constraints 
on the fact that we do not see a large population of H I clouds outside the disk of NGC 2997. The 
anomalous H I clouds within the disk, for which we calculate masses, become more consistent with 
those observed around M 31 and more consistent with the mass estimate for Complex C in the 
Milky Way. Despite this closer distance estimate, the difference between it and previous estimates 
is not enough to change the formation hypotheses for the individual clouds we identify. 



5. Conclusions 

Our observations of NGC 2297 reveal a galaxy with a l agging, thick disk and anamolous H I 



clouds in the halo simil ar to other galaxies, such as NGC 2403 (jFraternali et al.ll200ll ) and NGC 6946 



(jBoomsma et al.ll2008l ). The presence of a lagging disk, which may be caused by star formation 
activity, and the presence of large H I clouds and counter-rotating gas from accretion show that 
many processes contribute to the distribution of this cold gas. 

We propose that the lagging disk is the result of a galactic fountain process in which gas is 
blown out of the disk by supernovae. This gas experiences a decreasing gravitational force as it 
moves away from the mid-plane, so it drifts radially outward. As a result, this gas rotates in the 
same sense as the thin disk, but at slower velocities to conserve angular momentum. From IRAS 
60 and 100 /zm, and 1.4 GHz radio continuum fluxes we estimate fairly vigorous star formation 
rate in NGC 2997 of about 5 Mq yr^^ to support the renewal of this lagging disk. Of the total H I 
mass of the galaxy, 16 — 17% resides in this lagging disk. 

The most likely origin for the H I clouds we discover is accretion from a hot galactic halo or from 
nearby dark matter halos. The H I cloud masses are of order 10^ Mq. From their velocity offset 
with respect to nearby disk gas, we estimate a kinetic energy for the clouds relative to the rotation 
of the disk of ~ 10^^ ergs and conclude that it is unlikely their presence is due to supernovae. 
Meanwhile, the spur and tail may be the remnants of tidally disrupted gas rich dwarf galaxies. 
The anomalous H I we identify makes up approximately 1 — 2% of the total H I mass. Finally, 
the asymmetry of NGC 2997 to the northwest appears best explained by a tidal interaction with 
another galaxy, although it is not obvious whom the culprit may be. 

The features we find are not out of line with those observed in the Local Group, nor compared 
to other galaxies for which similar observations and modeling have been done, although the clouds of 
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NGC 2997 tend to be more massive than the clouds for which we have distance brackets in the Milky 
Way, even with an improved distance from the 2MASS Tully-Fisher relations that place NGC 2997 
12 — 23% closer. In spite of NGC 2997's relative isolation, we believe that a number of dynamic 
processes contribute to the sum of extraplanar H I gas, including star formation, accretion, and 
galaxy interactions. Deeper observations of NGC 2997 as have been accomplished with NGC 891 
would likely continue to reveal extended structure. 
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Fig. 1. — Channel maps of NGC 2997 covering 930.3 to 1029.3 kms^^ heliocentric velocity. Data 
is plotted at full spatial and spectral resolution. Contours are plotted for 3, 10, 20, 40, 60a, where 
(T = 1.9 X 10^^ cm~^ beam""^ channel"^. Grey scale is in units of mJy beam~^. 
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Fig. 4.— Top Left: DSS image of NGC 2997 ovcrlayed with H I contours. Top Right: H I total 
intensity map; contours are at 5, 20, 50, 80, 120a, where a = 1.9 x 10^^ cm~^ beam~^ channel"^. 
Grey scale is in units of Jy beam"^. Bottom Left: Intensity weighted velocity map of NGC 2997; 
contours range over 975-1200 kms"^ and are spaced every 25 km s ^. Bottom Right: the intensity 
weighted velocity dispersion map; contours are at 4, 10, 15, 20, 25, 40 kms~^. 
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Fig. 5. — The integrated H I profile of NGC 2997 over the full H I extent of the galaxy. The profile 
is centered on the systemic velocity of the galaxy, 1089 kms~^. 
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Fig. 6. — Left: A sample spectrum towards a position on the eastern side of NGC 2997 showing 
the thin disk rotation at ~ 120 kms~^; the asymmetry of the H I profile, characteristic of the thick 
disk, at ~ 60 kms~^; and anomalous H I at ~ —25 kms~^. Right: A spectrum from the western 
side of the galaxy illustrating the thin (~ 110 kms^^) and thick disk (~ 50 kms~^) components. 
The spectra are centered on the systemic velocity of the galaxy, 1089 kms~^. 
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Fig. 7. — The top panel is our rotation curve derived from a tilted ring model. It peaks at 226 km s~ 
around 12 kpc. The two bottom panels demonstrate the quality of our assumption that NGC 2997 
is an ideal, flat disk. Wc keep all other parameters fixed but allow the inclination (middle) and 
position angle (bottom) to vary — the rotational velocity for each tilted ring is fixed to the value 
found in the top panel. The declining position angle with radius suggests we were unable to ft a 
self-consistent rotation curve. 
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Fig. 8. — Rotation curves for both the thin disk (sohd line) and the thick disk (dashed Une). The 
average Av between the rotation curves is about 60 kms~^. 
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Fig. 9. — Limiting examples of our lagging disk models. Upper left- the best scale height and 
Zcyi, but too small a velocity gradient: there is not enough emission close to the systemic velocity. 
Upper right- the best velocity gradient and Zcyi, but no scale height: the emission is more spread 
out in velocity space, and the peak is lower in intensity and closer to the systemic velocity than 
seen in the data. Lower left- our best fit model with scale height = 0.8 kpc, Zcyi = 0.6 kpc, and 
velocity gradient = 38 kms~^ kpc~^. Lower right- the data. Contours are at 3, 10, 20, 40, 60a, 
where cr = 1.9 x 10^^ cm"^ beam~^ channel"^. 
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Fig. 10. — A total intensity H I map of NGC 2997 illustrating where slices have been made to 
analyze the anomalous H I content (see Figure [TT]) . 
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Fig. 11. — Select slices through the smoothed H I cube illustrate the extraplanar H I (see Figure [TOl 
& Table E]). (A) Clouds 1, 5, and the H I tail. (B) Cloud 2 with H I hole, the H I tail, spur, and 
additional examples of faint, counter-rotating gas. (C) Cloud 3 and H I spur. (D) Cloud 4 with 
H I hole, and H I spur. The positive direction corresponds to the western side of the galaxy. 



